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December) indieates that the annual deepening and subseaucnt decay of the 


easonal thermocline is accompanied by many random fluctuations, fhe 


Ci 


vhysical processes causing the decay and the fluctuations in the mixed- 
layer depth mare examined and qualitatively evaluated, 

The Tukey spectrum analysis programmed for the CDC 1604 electronic 
digital computer was used in analysis of oscillations at the botton of 
the mixed layer to determine the distribution of wave onergy for series 
of data taken at hourly intervals. An cnergy veal centered near 12 hours 
Was observed to predominate, and this energy is equivalent to taeG ot 2 
12-hour internal wave with height of 27.0 feet. The computer was also 
used to select the significant meteorological and oceanographic parameters 
which could be usec to predict a change in the mixed-layer depth by use 
of the BIND 07 multiple regression program. The most significant varameter 
was sea surface temperature. The best correlation coefficient for this 
parameter occurred for lags of zero to 12 hours. 

The results of the analysis support the theory that convection is 
the physical process which causes the seasonal decay of the thermo= 
cline curing the cooling season ani that short-term fluctuations of the 
mixed-layer depth are due primerily to internal waves. 

The authors wish to express their appreciation for the assistance 
and encouragement given then by Associate Professors G. H. Jung and 
J. B. Wickham, and by Professor D. A, Williams and the staff of the con 
puter center at the U. S. Naval Postgraduate School in this investigation. 
The authors also wish to thank the Pacific Oceanorraphic Groun of Canada 


for making the oce@mosmaphic data available. 
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Te troduction 

The vertical semperature structure of the near-surface layer in the 
open ocean is a complicated phenomenon which varies in spnrce and time, 
Certain aspects of the temperature structure are of great significance 
in determining sonar effectiveness. Thus, both the military and the com 
mercial fishing industries are interested in prediction of the thermal 
structure. 

Various processes continually modify the thermal structure to produce 
an essentially homogeneous (mixed) surface layer, underlain by a thin 
layer of rapidly increasing density, which varies in thickness through- 

ut the year. The nixedulayer depth (referred to as the MLD) is the depth 
below the water surface to which mixing has established isothermal condi~ 
tions, The lower boundary of the MLD is the top of the thermocline, a 

thin layer of large, negative vertical temperature gradicnt, usually associa~ 
ted with the layer of increasing density. <A complex intcraction of meteoro- 
logical, oceanographic, and (indirectly) astronomical factors provide the 
driving forewe for changes in the ocean thermal structure, 

Random and periodic fluctuations of the MLD have been examined by 
many investigators. It is known, as a result, that the MLD may vary fron 
a few inches in the heating season to several hundred feet in the cooling 
season. At times it may fluctuate tens of fect in a matter of minutes. 

The objective of this research is to identify more clearly significant 
processes Which affect the MLD during the cooling season and to relate the 
MLD to meteorological and oceanographic parameters, The processes capable 
of affecting the vertical temperature distribution are discussed.in the 
appendix. Certain of these processes will be Soneideee in this paper: 

(1) thermohaline poirettion (2) dynamic convection due to wind and wave 


action, and (3) unidentified processes which gencrate internal waves. 





Up to now, unly yualetytive descristions of the whxed-lnyor behavior 
during the cooling season are*iound in the 1 cerdvure, . Tiese usually 
state that increased wind stirrinc and surface cooling act tomether during 
thé cooling season to drive the seasonal thermocline down to its maximun 
annual depth. Our findings will show that, at Station Papa, thermohaline 
convection is the most important factor and that wind mixing (dynamic con- 
vection) is actually insignificant during the cooling season. 

Meteorological and oceanographic parameters contributing significantly 
to the change in the MLD during the cooling scason are determined, and from 
these the physical processes of importance are inferred. The investifa- 
tion was carried out with the ald of statistical correlation and weve. 
analysis programs, the BIMD 07 [16] and the Tukey spectrun analysis [15], 
utilizing the CDC 1604 digital computer. Regression equations and correla- 
tion coefficients are obtained by the BIND 07 program. One part of the 
Tukey spectrum analysis was used to determine cross spectra of the MLD 
fluctuations and the sea-surface temneratures, This vrogram was also used 

t 
taechne analysis of internal waves: 

The investigation is based on oceanographic and weather data obtained 
by the Fisheries Research Board of Canada, Pacific Oceanographic Group, 
for the months of October, November, and December during the years 1956, 


1957 and 1958 [6. 7]. 





Aa Backpround 
The attempt to identify the factors which bring about the hourly 
fluctuations in the MLD was based on the authors! earlier endeavor to 
derive an emirical method to predict the average }7.D on a daily basis 
during the cooling season. It became apparent then that the hourly fluc- 
tuations of the |iLD could be of equal or larger mapnitude than the avcragé 
deepening of the lL) during one or more days; thus, to use one or an 
average of two or more observed LD values per day as representative of 
the depth for any particular hour of that day is misleading eel The 
probable cnuse of these fluctuations is internal waves, 
Students at the U. S. Naval Postgraduate School during 1961 and 1962 
derived ompirical methods to predict the HI.D at Ocean Station Papa during 
ach of the 12 months of the year. In addition, the onset of the seasonal 
thermocline at Ocean Station Papa, which takes place in April or May, was 
studied by Clark in 1961 [2]. Clerk successfully releted tre thermocline 
nset to upper-air parameters. During the heating serson, when the sea- 
sonal thermocline gradually increases in magnitude, wind mixing processes 
have been clearly established as the dominant factor in determining the 
MLD, by Tabata (22) and Geary fio] in 1961. The cooling season remains to 


me scudicd. We have mde this attempt. 
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3. Deséra tion. of Ocean Otaticne aps 

The geographic lecation of Ocean Station Papa is 50N, 145W (figure 1). 
It is approximately 400 miles north of the boundary between the castern 
subarctic water mass and the subtropic water mass [sé] - the water depth 
at Papa's location is 2000 fathoms: thereforeg the thermal structure is 
in no way influenced by bottom topography. Also, its location is approxi- 
mately in thé center of thé Alaskan Current Gyral ke 9 ] (figure De 
Beeause of this location, the effects of advection are minimized and will 
not be considered. 

The weather ship normally maintains an on-station position within 
a ten-mile square centered on 50N, 145W. In order to avoid rough weather, 
the ship may stean weal eras within the 210-mile square grid (figure 2). 
Hach bathythernograph (BT) sounmling is identified by 2 two-lotter codec. 
Lt ce ship is on-station, the letter group OS indicates the position. 
The ship position for the greater part of the period of our concern wes 
Within a 30-mile square centcréed at OS. <Any dat& outside of this 30=mnil 


Square Were not consideréd in this study. 





f ovine op Tomer eae. Ah -Slanity SUG curtain. i Or Cool aie eae 
ficeu'e 3 reyrecents the dsea; of thes thsrrocline during the year 
1956. In August, tie curface Ilnyers have been heated to a maximum, approxi- 
mately 57F. Due primarily to wind ixime, there is a shallow mixed layer 
approximitely 100 feet in depth with a sharp underlying thermocline, With- 
in the thermocline, the water temperature drons to 4OF at a depth of 200 
feet. Progressing into the fall and winter, the MLD continually increases 
and the magnitude of the thermocline decreases. By the ond of December, 
bie nL) is neerly 330 feet in depth and the see-surface, tempcrature has 
decreased to approximately 42F,. Due to the presence of an intense halo- 
elaine, tne MLD reaches a limiting depth of 330 feet in Januarytand 1s0— 
thermal conditions prevail down to this depth for the remainder of the 
Wanwer . 

Thewbasic salinity structure for the eastern subarctic Pacific water 
mass is shown in figure 4. Three zones’ are indicated. The upper zone 
extends from the sea surface to a depth of 330 fect, and is characterized 
by relatively lowesalinity water (32. Tye [21] - Pies aad 


° 


Zone 1S isonaline and is mar%ed by the presence of ths sonsunal thormo— 


cline. JUecause of the isohaline conditions in the unner 330 feet, we 


e 


@scume thal sicnificant variation in density is governed by changes in 
temperature only. 

the halocline represents a transition zono between the uyper and 
dower zones. here, the salinity increases markedly with depth (as much 
as one part per thousand within an interv2l of 330 feet), The halocline 
Represents a limiting depth Lor the seasonal thermocline. Tt is si:ntlar 
to having an ocean bottom at this depth as far as convective mixing is 


concérned. The density incréase®is so great throurch this layer that missing 


cannot take place through it [20]. Within the lower zone (> 650 fest) 
ne: salami magautlivy inereices wet -aepe: tO the coean bobcon, 
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oy heen 

The purpose of this statistical study of the MLD is three=fold, 

An attempt is made: (1) to determine the relative importanec of thermo- 
haline conveetive mixing and wind-induced mixing on the deepening of the 
thermocline during the cooling season; (2) to determine the lag between 

a signifieant meteorological or oceanographic event and a subsequent change 
in the LD; and (3) to examine the importanee of the seemingly randon 
snorteterm fluctuations of the MLD and to provide a description of those 
internal waves present. 

First, the parameter or parameters are identified which have the 
largest linear correlation with the deepening of the MLD. Next, a further 
investigation has been undertaken of these significant parameters to deter-= 
mine the lag in the mixing process, that is, the time between the onset 
of a change in the parameter and a subsequent deepening of the MLD. 
Finally, the observed variations in the MLD which eannot be accounted for 


by either mixing process are investigated, 


5,6. Method of Investigation 

The parameters to be investigated are sub-divided into three groups, 
Those related to convective stirring are: (1) temperature of the sea, 
at the time of the sounding and with a 2-hour lag; (2) air temperature, 
at the time of the sounding and with a 2-hour lag; (3) devwnoint tempera- 
ture and wet-bulb temperature, at the time of the sounding; (4) relative 
hunidity, at 2-hour end 5~hour lags; and (5) low-cloud cover, with a 2- 
hour lag. 

Parameters chosen to represent wind and wave mixing are: (1) wind 
speed, at 2-hour and 5—hour lags; and (2) height of significant waves, 


at 2-hour and 5-hour lags. 


6; 





Parameters chosen which were shown to have 1 correlation with the 
thormal structure by others in previous studics, but did not conveniently 
fall into the above categories, are: (1) 500-mb height, in geopotential 
meters; (2) 500-mb temperature and 500-mb relative hwnidity, at a 2-hour 
lag; and (3) surface-pressure change, with a 2~hour and 5-hour lag. 

These data were obtained from synoptic 3—hourly reports, twice-daily 
constantepressure records, and the BT traces. The available BT's for the 
period of investigation were tabulated in two groups, 02004 and 17004, 
for cach year to climinate any possible yearly anomalies or diurnal trends. 
An investigation by Scripps Institution of Oceanography of the annual MLD 
fluctuation in the northeastern Pacific Ocean shows that the MLD most fre- 
quently ohserved in a given month often varies considerably from yeer to 
year [16]. The available literature on diurnal trends during the cooling 
season is not so clear. A preliminary investigation of the data showed 
no marked diurnal trend. This indication was supported by Dr. Tully of 
the Pacific Oceanographic Group, in a conversation with the authors, in 
which he stated thai during the cooling season the MLD is toc deep to have 
a clearly defined diurnal fluctuation, 

With the aid of a multiplesregression and correlation analysis, the 
BIMD O07 programmed for the CDC 1604 high-speed computer, the linear- 
correlation coefficients were obtained for the pre-selected group of me- 
teorological and oceanographic parameters with the MLD. The BIND 07 pro- 
gram was used because it can select different sub<-samples of data that 
are obtained from the same population. Thus, the program would perform 
linenr-correlation analysis and linear multiple regression on the data 
for 02004 and 17008, separately and in coribination. In this way, a close 
examination of possible diurnal trends could be made for each of the three 


years studied, 1956 through 1958, The analysis provided regression equations 





for nine different data groups for the periods shown in Table I, 

The BIMD 07 program yields best results when the number of indepen- 
dent variables is not greater than one-half the number of values of the 
dependent variable. Since the smallest number of values for the depen- 
dent variable is 36 (at 02002, 1958)!, only 18 independent variables were 
investigated. A set of 18 secondary parameters was tabulated in the event 
that the analysis of the 18 primary paraneters proved inconclusive, 

There was no attempt by the authors to eliminate possible duplication 
of parameters that are indicative of the same process, For example, dew 
point temperature, wet—bulb temperature, relative humidity, and low-cloud 
cover are all related to evaporation; but it was left to the BIMD 07 program 
to discern which-ones, if any, were significant in contributing to the 
prediction of the deepening of the IiLD, 

The BIMD 07 is unable to produce a miltiple linecar-regression equation 
if any independent variable is a linear combination of two or more inde- 
pendent variables. Therefore, temperature of the air minus temperature 
of the sea, which is often used as an indicator of thermohaline convection, 
could not be investigated in conjunction with the air temperature and sea 
temperature, Gonsequently, three additional problems were computed using 
the air-minusesee temperature parameter in place of sea=temperature para- 
meters. voce Table II for the tabulated results of thé*1l2 problems. 

Tntenalysis of all thé data, certain facts and tendencies are evident; 
these are discussed in the following two paragraphs. 

The only parameters that show a consistently high linear correlation 
with the 10.0 in every sample or sub-sample are the sea-surface tempera 


tures, at the time of  e sounding and with a 2-hour lag. They also accounted 


1(This small number results because BT data for October, 1958 are jusse 
ing; only a small number of observations was available in 1957 as well, since 
December, 195% observations were not taken. 


8 





for almost all of the variance explained by the regression equation using 
11118 varameters. Figure 4 shows the close relationship between the 
scea.surface temperature and the MLD from October through December, 1958. 
Likewise, Zottel [25] suggested that sea-surface temmerature might be a ? 
good "indic&tor" for thé LD. 

Wind speed and height of the significant waves, the parameters chosen 
as indicators of wind stirring, have either very low negative or low posi- 
tive correlation coefficients, 

In order to determine the time required for an increase in convective 
activity due to a reduction of sea-surface temperature to effect a sub- 
sequent deepening of the MLD, the MLD's for 17002, 1956, were tabulated 
with the sea temperature at lags of zero, two, five, eight, eleven, four- 
teen, seventecn, twenty-one, and twenty-three hours as the indevendent 
variables. A simple linear-correlation program was used to determine 
that the lag yielding the best correlation was two hours (Table III). 

Also, two one-hourly series of MLD for 79 hours in Oetober, 1957, and 
for 51 hours in November, 1958, were tabulated with their corresponding 
sea~surface tomperatures. The energy-density spectrum and the cross=spec- 
trum of both the MLD and the sea-surfnaec-temperature time series was ob- 
tained from the Tukey analysis. Maximum lags of 39 hours and 24 hours 
were used for tne 79-hour and 5l-hour series respectively. The two co- 
spectra for the 79-hour and the 5l-hour series show large values for 
h's around 6 hours and 4 hours respectively (figure 5), 

From the results of the simple linear ztzorrelation and the Tukey analysis,- 
it is seen that all lags from two tosix hours show high correlation. Re- 
yond 12 hours the correlation rapidly decreases, 

A linsarerepression equation was derived for 17004, 1956, using sea 


tomperature with a 2~hour lag as the single independent variable (Table IV ), 


1 on eross-correlation, not shown 
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The fact that the results of this equation as measured by the standard 
error comoare payee ot with the multipleregression equation derived from 
13 independent variables for the same year supports the initial analysis 
of the problem. If a regression equation were to be used to forecast MLD, 
one using only the sca-surface temperature taken with a lag of up to l2 
hours would be much more reliable than an equation with more parameters. 
This is true simply because it would have a larger number of degrees of 
freedom. Note that all of the regression equations have a standard orror 
of from 15 to 30 feet, with an aversge for all problems of 22.15 feet 
(Table V, Table VI), 

There are several possible explanations for this standard error of 
more than 20 feet: (1) it is possible that not enough parameters were 
used, and that by substituting for primary parameters which showed a very 
poor correlation, some of the secondary parameters or primary parameters 
with different lags, the standard error of the new regression equation 
would be reduced; or (2) there is a strong possibility that a non-linea 
relationship between a parameter and the MLD could account for a large 
part of the standard error (however, the BIMD 07 program cannot detect 
this); or (3) the error may be due to internal waves. Of these possi- 
bilities, (3) seems the most likely reason to explain tho observed stan- 


dard crror, 


LO 





( 


Ox Analysis of Internal Waves | 

The base of the ji2iM@d layer in the eccan 1s subjcct to nany vertzcal 
oscillations. Figure 7 shows a series of daily MLD's for the period Octo- 
ber through December, 1956. ach vlotted point on the curve is an average 
daily value of the 11.D, Two individual soundings are normally used; how- 
ever, un to 24 soundings were available on several days. It is observed 
that the MLD undergoes a seasonal downward trend. During the October 
through December period, the seasonal trend of the MLD finds it increas- 
ine from approximately 130 feet to 330 feet. This is an av@rage detpen-— 
ing of 2.17 fect per day. Also, barge vertical fluctuations ames super= 
imposed on this downward trend, which become larger in magnitude as the 
mixed-layer thickness increases. Note in figure 7 that the vertical fluc- 
tuations of the MLD can be ten times greater than the average daily deepen- 
ing. In addition, it is seen (figure $) that hourly fluctuations can be 
of equal or greater magnitude than the average daily deepening during the 
al. 

Fisure & represents a one-hourly time series for 51 consecutive hours 
during November, 1958 (series 1). A second one~hourly time series rep- 
resents 79 consecutive observations taken during October, 1957 (series 2). 
The magnitude of the MLD fluctuations is about the same, whether it ic 
computed from hourly, twice-daily, or daily obsérvations. The ®avéragt 
fluctuation per hour of the MLD for series 1 is 18 feet. The average 
fluctuation between the consecutive 02004 and 17004 soundings for the month 
of November, 1958, is 21 feet. Also, the averaged daily fluctuations during 
November is 25 feet. The numbers for series 2 are sinilar. A December, 
1956, one~hourly time series (figure 9) was available but was too short 
for detailed analysis. However, it illustrates the extreme fluctuation 


of the MLD which takes place over a short time interval. Over a 3-—hour 


ll 
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orilody tial ILD incr+3sed. from 2/0 fea to 5/0 Toot, anc unen Uceres ca 
Ale LOnc J) Cet. 

Since mixing »nrocesses cin account for only an increase in the ILD, 
and since observed fluctuations are in both directions (up and down), 
and are a consistent feature of the data, these observed fluctuations 
are considered to be due to internal waves. 

Internal waves may occur within stratified water and in water in 
which the density increases with depth. The largest vertical displace- 
ments of the water particles are to be found within the boundary between 
layers of different density, The displacement amplitude diminishes above 
and below this boundary, approaching zero at the sea surface and the ocean 
bottom. The density difference between the waters above and below the 
seasonal thermocline is, of course, much less than the density difference 
between the air and the water at the sea surface. Therefore, a boundary 
Within the ocean can be much more easily displaced than the surface of 
the sea [17] Ufford [24] states that surface waves require 30,000 tines 
more energy to start and maintain than is required for internal waves of 
the same amplitude, Lnergy for internal waves could come from surface 
disturbaneés or from currénhts within the oc@éan, : 

Besides rapid vertical oscillations in the MLD, it was observed that 
the thickness (or vertical extent) of the thermocline layer below the 
MLD may vary markedly with time. Fisure 10 shows the thermocline thick- 
ness for scries 7. Variations in the thermocline thickness indicate that 
intvernal-wave characteristics of phase and emplitud@ very with depth. 
Also, the temperature gradient in the thermocline varies markedly because 
of changes in its thickness. The actual terporeturc difference through 
tne tnermocline is 2 roletivel:: stable Beature and is therefore predic= 
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Wackucsa to skady-zepeach time cerios. The magngesiccumeas to [inc guy 
Orecovacint geriods ane 2scocizmed amolitudcs recent) in tice oscil 
dations. The, Tukey anilysis estimates the cnucorcyedensit; as @ function 
of porind in a siven series. This is done by commnting the wto-correla~ 
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Spoeurn. ectir-te is obtained fron the Fourier cocint transform of eine 
nuto—corrélation function, 

Initially, 2 correlorram vas plotted for both time series (fisure 11), 
ie corrélorranm iss@ plot of the non-normalized auto-correlation coerri-~ 
cients as a function of the lag (k}. Inspection of this correlogzran appears 
wevindieete that the svectrun of series @ has emelletively large enerpy- 


nd a cecondary maximum with a 


oF 
fo 


acne2 ty maximum meer the 12-heur @crio 
peak near six hors; and that the spectrum of series 2 has only one maximun 
Rear 12 hours. 

Next, the energy—density spectrum of both time series was obtained 
from the Tukey analysis. Figure 12is a plot of the energy—density spectrun 
for series 1. A maximum lag (mn) of 24 hours was used, which allowed the 
spectrun to includ@ a range of periods from two hours through 48 hours. ~- 
The snectrim helns confirm the indications of the correlogram, A maximum 
of energy of the variance of the M.D is shown to be produced by fluctua- 


tions with periods on the order of 12 hours. Also, a large energy—density 


is indicated for periods arcund six hours, There is an energy gap for the 


ike 





periods around eight hours, and the cnergy for scriods less than five hours 
areops off rapidly to@@n insignificant amount, 

The totil energy is a known quantity, cqual to twice the variance, 
By measuring total area under the spectrum curve and the area under each 
frequency peak between half-energy points, equivalent amplitudes for waves 
with neriods near 12 hours and six hours were Sceeel, These computed 
amplitudes are 18 feet and 12 feet for the 12-hour and the 6~hour periods, 
Gespectivey. 

The energy-density spectrum for series 2 is shown in figure 13, A 
naximun lag (m) of 30 hours was used, which allowed inspection of periods 
ranging from two hours through 60 hours. A relatively large cnergy maxi- 
mum is again indicated for periods arcund 12 hours. ‘The curve then falls 
off rapidly, but shows a slight maximum near the 6-hour range of periods. 
A single sine wave with energy equivalent to that under the 12-hour peak 
would have an amplitude of nine feet. 

This analysis reveals three features of the thermal structure at 
Ocean Station Papa during the time interval investigated. 

(1) Yortical oscillations of the MLD are present at the thermocline 
during the cooling season, ‘Tne magnitudes of these oscillations are 
similar in all data, whether they be from hourly, twice-daily, or daily 
observations. 

(2) Although the time series available was mich too short on which 
to base strong conclusions, the analysis did reveal that the mijor energy 
peaks occurred with a period of approximately 12 hours with an average 
amplitude of 13.5 feet. 

(3) The value of the standard error of estimate resulting from the 


BILD O7 regression equations is very close to the height of the postulated 


12-hour wave (that is, 27 feet). Therefore, it is concluded that there 
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is a direct relation botween this error and 


waves durincs the period of investigation. 
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as Conclusions and Recomiendations 

As a result of this stud; of the influences on the MLD during the 
cooling season, the following conclusions can be reachea. 

(1) Thermohaline convection is the important process for the deepen- 
ing of the MLD in the autumn season. 

(2) During the period of study, the MLD is too deep to be affected 
by Wave stirring or diurnal trénds. 

(3) The time lag between surface cooling and deepening of the MLD 
by the convective process is small; this suggests that once the cooling 
season begins, convective circulation is essentially continuous, This 
circulation may be maintained during periods of near-zero heat flux across 
Le, eir-sea interface by wind energy, as suggested by Jakond [23] ; hove 
ever, until an injection of notential energy in the form of rapid surface 
ecoling, the MLD will remain relatively static. If the increased surface 
density caused by surface cooling is sufficient, the convective circula- 
tion intensifies and drives the MLD deeper. 

(4) A major portion of the MLD variability appears to be associated 
with internal waves. The analysis of internal waves revealed that the 
major enerry peaks occurred with a period of approximatély 12 hours with 
en average amplitude of 13.5 feet, 

Barly in 1953, at a Navy-sponsored thermocline conference attended 
by 20 of this country's most prominent oceanographers and meteorologists, 
the following question was asked his colleagues by Dr. T. F. Malone: "How 
accurate a prediction of the I1iLD and its hourly fluctuations is required 
by the nilitary?" [19] This remains a vital qué@stion. 

Tf the military is willing to accent forecasts of the MLD within 
ranges of 20 or 30 feet, then existing methods of prediction will suffice. 


Any increase in accuracy over cxisting schemes for predicting the MLD at 
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1 particulier hour on a narticular dey can be accomolished in areas whene 
internal waves are important only after prediction of the internal wave 
pattern is possible. This is a complicated problen, 

Ifa sufficiently lonz, detailed timé Series of BT's could be made 
evailable, internal waves could be filteréd out by use of a technique such 
as proposed by Linnette EAR This would leave a time serics of the MLD 
that is dué entirely to convéctive mixing or other causes, thus enabling 
a mich better quantitative analysis of the convective process, This is 
suggested as the next logical extension of the research reported in this 


paper. 
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Table I 


1956 1957 1958 
02008 = 48 es 26 
17008 61 16 13 

109 90 19 


(0200% is 1600 local mean time) 
(1700K is 0700 local mean time) 


Bathythermograph Soundings Available 
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Table IV 


Linearecorrelation cocfficient 
(ior ioe) leo 


Standard error of the estimte 30.33 


Range of residuals 2 


Regression equation: MLDesz, = 1152.94 = 20.24 (ia) 


Regression Results for One Independent Variable 
(17004, 1956) 
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Multiple-Regression Results for 18 Independent Variables 
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Table VI 


Multiple-correlation coefficient 
(for 18 independent variables) «596 


Linear-correclation coefficicrnt 


(Por Tevand Tees) -.ol4 
Standard error of the estimate Zoos 
Range of residuals tea 


Standard deviation L350 


Three-Year Averages for Regression Results 
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Heat Transport in the Sea and Exchange across the Air-Scea Interface 
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A. Thermohaline convection. The vertical displacement of srall par- 
cels of water occurs when a small vart of @ Water mass is neaviecr than tae 
weter underneath it. To restore this unstable Gondition to equilibrium, 
the heavier water tends to sink while the lighter water rises. Associated 
with these forced vertical movements of small parcels there is also a trans-— 
port of the characteristic properties of sea water in a vertical direction 
which leads to an equalizetion of any vertical differences in these proper- 
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ties which may be present. It is believed that convective mixing 
important an agent of heat transfer in the occan as wind stirring during 
the cooling season (2. Lay 20] : 

An initial increase in the density of small particles at the surface 
accompanies an increase in salinity (due to evaporation, or to the form- 
tion of ice) or a decrease in temperature. One or more of these in combi- 
nation may be involved in thermohaline convection. In lower latitudes, 
where there are anly small variations in the temperature, the heat loss 
1s outweighed by the effect of evaporation; in polar regions, in addition 
to radiation and evaporation the increase in salinity due to formation of 
ice is also effective. In temperate latitudes the heat loss by radiation 
temivaice decisive: fmctor. 

Convective sensible heat exchange (conduction) between air and sea 
surface has been studied by Kuhbrodt [4] for the North Atlantic. He has 
shown the convectional heat flux to be approximately 20 on cal/en per day 
during the period of the year when the air is cooler then the sea surface, 
This value, even as a rough estimate, is too large to be neglected in the 
heat budget of the ocean. If the surface of the sea is warmer than the 


atmosphere, the air is heated at the interface and the vertical stratifi- 


38 








7 


cation of the air becomes unstable; the air becories turbulent with a con 
stant replacement by cool air at the sea surface and the vertical heat 
transport becomes large. 

Values given by Defant [4] for the meet, Loss 2m tints. of pr sal/en@ 
per day for 50N latitude are 116 units due to effective back radiation, 

78 units due to evaporation, and 20 units due to convection of sensible 
heat (conduction) as a yearly average. Gertainly, during the cooling 
season heat loss due to conduction would be even more significant. These 
values are derived, assuming that the heat exchange through the ocean sur- 
face occurs independently for each separate latitude belt. Therefore, no 
neridional heat exchange (by ocean currents and by horizontal mixing) was 
considered. 

B. Dynamic convection (forced vertical mixing due to wind and wave 
action). This process has been well studied and considered to be the 
dominant process of vertical mixing in the heating season. Laevastu 12] 
has proposed an empirical formula for determining the MLD by significant 
wave height alone. 

C,; Dynmemie convection due to forced vertical mixing of ogean currents 
can result in convective mixing, and influences the ML, This 
process was not considered because of the difficulty of measurenent; hovw- 
ever, Since the Station is located in the center of the Alaskan Gyral, 
ecean-current effects are vrobably small. 

D. Thermel conductivity can take place when a vertical temperature 
Gradient exists in the ocear, Heat is transferred by molecular heat con 
Guction processes. Defant [4] eoneludes that this process is insignificant 
in oceanozranhic investigations due to the extremely long period involved 
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in the conduction process. Therefore, this process will not be considered 
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7 wnbtrainnent due to turbulence at tue Cycnecline. Durin: the 








coolinjss scason in the open ocean, a well-defined isothormal surface layer 
esiecs Whee extends downward. to ane teri 1cc ole rapid (ecresa2tn si icr 
perature. Investizations of the subarctic water indicate an isohaline 
condition to a denth of 23) feet, so the cffect of salinity on density 

can be ignored (23h and thus we nay consider the thermocline as a2 pycno= 
cline. Because less dense warner water overlics dense cold water, we have 
a stable condition mich like a surface inversion in the atmosnhere. When 
surface cooling takes place, thereby creating an unstable condition, free 
eonvection ean) occur. Cromvell[3] has demonstrated that this free convec- 
tion in a mocel causes a turbulent exchange across the thernocline which 
is one-way. tne fluid particles which move upward from the region below 
the thermocline are rapidly deformed and mixed throughout the upper leyers; 
fluid particles which move downward toward the quiet layer below the MLD 
are buoyed vpward intact. Thus, the upper layer increases in thickness, 
but decreases in temperature,it the expense of the lower layer. This same 
phenomenon is observed when the turbulence is induced by wave action or in- 


bermal waive action [2] * 
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